hep-th/0202067 
KCL-MTH-02-03 



Boundary states for WZW models 



Matthias R. Gaberdiel* 

Department of Mathematics, King's College London 
Strand, London WC2R 2LS, U.K. 



and 



Terry Gannon* 

Department of Mathematical Sciences, University of Alberta 
Edmonton, Alberta, Canada, T6G 2G1 



Abstract 

The boundary states for a certain class of WZW models are determined. The models 
include all modular invariants that are associated to a symmetry of the unextended Dynkin 
diagram. Explicit formulae for the boundary state coefficients are given in each case, and 
a number of properties of the corresponding NIM-reps are derived. 



February, 2002 



* mrg@mth.kcl.ac.uk 

* tgarmori@math . ualberta . ca 



1. Introduction 

Given a conformal field theory defined on closed Riemann surfaces (a closed string 
theory) we can ask whether the theory can be defined on Riemann surfaces with boundaries. 
More precisely we can ask which boundary conditions can be imposed at the various 
boundaries. From a string theory point of view, this is the question of which open strings 
can be consistently added to a given closed string theory. In general, rather little is known 
about when a 'complete' set of boundary conditions can be found. It is believed that 
for closed string theories that are consistently defined on arbitrary genus surfaces such a 
solution exists, but not even this is clear. Even if the higher genus property is sufficient, 
it seems unlikely to be necessary. 

One of the conditions that every consistent set of boundary states has to satisfy is the 
so-called Cardy condition [0 : suppose that the conformal field theory is a rational confor- 
mal field theory with respect to some symmetry algebra A, which possesses iV irreducible 
highest weight representations. The Cardy condition then implies that every set of M 
boundary conditions determines N M x M matrices of non-negative integers, one for each 
highest weight representation of A. Furthermore, if the set of boundary conditions is com- 
plete in a suitable sense ||, these matrices actually form a (Non-negative Integer Matrix) 
representation, or NIM-rep for short, of the fusion algebra (see f^fill f° r an introduction 
to these matters). 

The condition of a representation of the fusion algebra to be a NIM-rep is quite 
restrictive, and one can therefore attempt to classify all such representations (irrespective 
of whether they arise from boundary states or not). If modular invariance is indeed a 
sufficient condition for the construction of the boundary states, the classification of NIM- 
reps then gives a restriction on the possible modular invariant theories. This program has 
been performed for the case of the WZW models corresponding to su(2) in ||. However, 
in general this approach does not seem to be powerful since there are typically many more 
NIM-reps than modular invariant conformal field theories ||, the best known examples 
being the tadpoles of su(2) (There are also many modular invariant partition functions 
that do not have a NIM-rep H; however, at least some of them do not define consistent 
conformal field theories, and it is therefore conceivable that modular invariance is indeed 
sufficient for the construction of a boundary theory.) In the following we shall call a NIM- 
rep 'physical' if it is compatible with a modular invariant (see the end of section 2 below) . 
If the NIM-rep possesses in addition compatible structures corresponding to the defect 
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lines and 3j- and 6j-symbols (see for example 0), we shall call it 'fully realised'. Fully 
realised NIM-reps can be obtained from a subfactor . 

It is generally believed that the conformal field theory whose space of states is de- 
scribed by the diagonal modular invariant defines a consistent conformal field theory that is 
defined on arbitrary Riemann surfaces.* Given this theory, we can obtain other consistent 
modular invariant conformal field theories by applying a (global) symmetry transforma- 
tion to the left-moving degrees of freedom, say. For example, for the case of the WZW 
models corresponding to the simple Lie algebra g, any symmetry of the (unextended) 
Dynkin diagram can be employed. Symmetries of the unextended Dynkin diagram exist 
for g = su(n) and Eq, where they correspond to 'charge conjugation'; for g = so(2n), where 
they correspond to the 'chirality flip'; and for g — so(8) where there is in addition 'triality'. 

The boundary states of the usual diagonal modular invariant theory are known in 
detail and explicitly (they are described by Cardy's formula [Q), but the boundary states 
of the theories that can be obtained using one of the above symmetries are less well 
understood. In this paper we shall make a proposal for the boundary states of these 
theories. Our construction is based on the observation that the problem of constructing 
boundary states for these theories is essentially equivalent to the problem of finding the 
twisted boundary states of the original theory. These twisted boundary states satisfy two 
sets of Cardy-like conditions, one coming from the overlaps between the twisted boundary 
states among themselves, and one from the overlaps between twisted boundary states and 
the usual Cardy states. These two Cardy conditions combine to give a rather powerful 
constraint on the twisted boundary states, that will allow us to make a very suggestive 
proposal. We are also able to show that the NIM-rep entries it defines are always integers, 
and we have checked explicitly for numerous examples, that these integers are in fact 
non-negative. 

From a mathematical point of view, the WZW models are particularly interesting 
since many of its quantities can be given a natural algebraic interpretation in terms of the 
associated affine untwisted Kac-Moody algebra g. For instance, the highest weight repre- 
sentations of the symmetry algebra A are labelled by the integrable highest weights A G P+ 
of g. The fusion rules N% can be interpreted as Weyl-folded tensor product coefficients 

* In the diagonal modular invariant, left- and right-moving representations occur in conjugate 
pairs. See section 2 for a more detailed definition. 
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of the underlying finite-dimensional algebra g. And the matrix S which diagonalises the 
fusion rules describes the r i-> -1/r transformation of the characters X\( r ) °f fl- 
it is therefore natural to ask three questions about the NIM-reps in WZW models. Namely, 
find algebraic interpretations for 

(a) the boundary state labels x G £>, that is to say the rows and columns of the 
NIM-reps Af\ (in the type III sub factor language these are the M — N 
sectors) ; 

(b) the NIM-rep coefficients themselves; and 

(c) the matrix ip diagonalising the matrices Af\, and expressing the boundary 
states as linear combinations of the Ishibashi states. 

In this paper we suggest answers for the NIM-reps that arise as described above, and that 
are therefore associated to the symmetries u of the (unextended) Dynkin diagram of g. In 
particular, 

(a) the x G B are integrable highest weights for the twisted affine algebra g w , 
or equivalently the highest weights of w-twisted representations of g; 

(b) the Af^ x are the coefficients of the fusion of the w-twisted g-representation 
x with the untwisted g-representation A; and 

(c) ijj is the matrix S describing how the characters of the twisted algebra g w 
transform under r i— > —1/r. 

The exponents fx G £ can also be interpreted as the highest weights of another twisted 
algebra, the so-called orbit algebra g. Just as the symmetries of the (extended) Dynkin 
diagram of g give rise to simple currents and thus to symmetries and gradings of the fusion 
coefficients, we will find that the symmetries of the Dynkin diagrams of g^ and g give rise 
to symmetries and gradings of the NIM-rep coefficients. 

We also interpret our answers to (a), (b) and (c) in terms of data for certain affine 
untwisted algebras; in particular we can express our NIM-rep coefficients using ordinary 
fusion coefficients, proving the integrality of our NIM-reps. 

Answers to these questions, which appear to be similar to ours (although they were 
obtained by independent arguments), appear in work by Fuchs, Schweigert, and collabo- 
rators (see in particular ||,|10]]). However, we feel that the fundamental simplicity of the 
above algebraic descriptions, and the fact that one can be completely explicit, is rather 
elusive in their work. In particular, their description is in terms of an orbifold construction, 
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and this hides the simple underlying structure we exhibit in our paper. They have also 
not yet shown that their conjectured NIM-reps are in fact integral. 

Given this algebraic description of the NIM-reps, it seems plausible that our con- 
struction may generalise further to symmetries that are only symmetries of the extended 
Dynkin diagram (and that therefore give rise to simple current modular invariants). This 
is illustrated in section 6 where we construct the NIM-reps for the simple current modular 
invariants of su(p) with p > 2 prime (for p = 2 we ignore the easy case k/2 odd). As long 
as p 2 does not divide k + p we can furthermore prove that the NIM-rep entries are indeed 
non-negative integers. A more complete description of these 'simple current' NIM-reps and 
their non-negativity and integrality will be given elsewhere. These simple current NIM- 
reps have been addressed in [JTJ in considerable generality, although in our opinion their 
description is much less explicit than our section 6. They also provide an argument for the 
integrality and non-negativity of the NIM-rep coefficients, but it depends on the validity 
of their conjectured formula for the fusion coefficients of the chiral algebra extended by 
simple currents, and this conjecture has not yet to our knowledge been rigorously estab- 
lished. For recent work on a completely different approach to proving this integrality and 



non-negativity, see 12 



Surprisingly little is known about explicit NIM-reps, even for the WZW models. The 
known NIM-reps include at least one corresponding to each modular invariant of su(2) [|13| 
and su(3) [|13],|7],|^,|||3| — see @ for the lists; at least one NIM-rep for each known modular 
invariant of su(4) |l|,0,|J|,|n|; and a few NIM-reps for su(5) and G 2 |]|0,[^g,|IlJ . The 
list of NIM-reps for su(2) is complete |T3[]; the lists of 'fully realised' NIM-reps for su(3) 
and su(4) have also been claimed to be complete [§ , although the arguments have not yet 
appeared in print. On the other hand, it is known that the lists of all NIM-reps (without 
assuming them to be 'physical' or 'fully realised') are incomplete for all algebras other 
than su(2). As mentioned before, NIM-reps associated to conjugation automorphisms and 
simple currents are constructed for instance in PJll]]; however their formulae are far more 
complicated and less explicit than what is done here and in other papers quoted above. 
In the type III subfactor framework, NIM-reps compatible with the modular invariants 
corresponding to simple currents |7] and conformal embeddings |T5|,[7[1, were shown to ex- 
ist for su(n), and thus these NIM-reps are in principle computable from the subfactor 
machinery |16|J7[|, although this computation is not entirely straightforward and has not 
been performed. The diagonalising matrix ijj, however, is at present inaccessible from the 
subfactor approach. 
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While this paper was being completed, [I7J appeared in which the NIM-reps corre- 



sponding to charge conjugation C for su(n) are obtained, using different methods than 
ours. (In particular, their analysis does not involve twisted affine algebras.) Their NIM- 
reps seem to be consistent with ours, although their diagonalising ^-matrices differ from 
ours by some (immaterial) signs, and their boundary state labels are different {e.g. for 
su(5) level 2, their boundary states (mi, 7712) = (1,1), (1,3) and (2,1) correspond to our 
labels (a±, (12) = (0, 1), (0, 0) and (1, 0), respectively). After completing this paper we were 
also made aware of the paper [ff8[ whose results have some overlap with our section 3. 



The paper is organised as follows. In section 2 we describe our conventions and some 
background material. Our construction of the NIM-reps is motivated in section 3 where 
we explain how it can be understood in terms of constructing twisted boundary conditions 
for the original diagonal theory. In section 4 the NIM-reps are explicitly described, and 
their relation to ordinary fusion matrices is exhibited, thereby proving the integrality of 
the NIM-rep coefficients. Details of the underlying calculations are given in section 5. In 
section 6 we describe the NIM-reps for some simple current modular invariants, and prove 
that, under suitable conditions, their entries are indeed non-negative integers. Finally, 
section 7 contains some conclusions. We have included an appendix in which the graphs 
of some of the NIM-reps are explicitly given for small rank and level. 



2. Conventions and background material 

We will be interested here in the rational conformal field theory associated to an affine 
untwisted Kac-Moody algebra g at positive integer level k. Let us denote by P + = -P+(g) 
the set of all (finitely many) integrable level k highest weight representations of q. The 
corresponding highest weights can be written as A = J27=o ^i-^i = (^o> Ai, . . . , A n ), where 
the non-negative integers A^ are the Dynkin labels of A, and the Aj are the fundamental 
weights of g. For instance, the vacuum representation corresponds to the weight kAg. 

The characters X\( T )i ^ £ P+, °f a rational conformal field theory carry a represen- 



tation of the modular group jnj, and in particular 



Xa(-Vt)= ^aX^(t), (2.1) 



where S is a unitary symmetric matrix. By the Verlinde formula [[2DJ the S'-matrix deter- 
mines the fusion rule coefficients of the theory as 



000* 

Nx » - ^ S 0K ■ (2 ' 2) 



Charge conjugation C is a permutation of -P+, sending the weight A to the weight A c 
contragredient to it; it is given by C = S 2 , and obeys iV° = 5 Mi ca- C always corresponds 
to a symmetry of the unextended Dynkin diagram (although this symmetry may be trivial). 

Another permutation of P+ which is of fundamental importance in the theory, is due 
to the simple currents (see for example [ PT|,P2| ). These are the weights J G P+ with 
quantum dimension = 1; their fusion rules Nj^ = d^j^ define a permutation \x i— > J \x. 
This permutation obeys 

Sxj» = e 27r[Q ^S Xfl (2.3) 

for a rational phase Q(A). For all WZW theories except Eg level 2, the simple currents 
correspond to symmetries of the extended Dynkin diagram. Simple currents give rise to 
symmetries and gradings of fusion coefficients 

^a M 7^0 ^Q(X) + Q(fi) = Q(u) (modi). 

We will encounter simple currents in our formulae in section 4 for the NIM-rep coefficients, 
and more explicitly in section 6. 

We shall only consider conformal field theories for which the left- and right-moving 
chiral algebra is the same, namely the untwisted Kac-Moody algebra q. The space of states 
of the full conformal field theory can then be decomposed into representations of two copies 
of g, and thus can be written as 

H= M Xfl Hx®H^, (2.5) 

where are non-negative integers that describe the multiplicity with which the various 
tensor products of representations appear in H. We shall call the theory defined by M = I 
the 'diagonal' theory; from the point of view of constructing boundary states, this is the 
simplest case (as we shall see momentarily). We shall mainly be interested in theories 
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that can be denned on the torus; this requires, in particular, that the partition function 
corresponding to (|2.5| ), 

Z(r)= Yl M vX a(t)x m (t)*, (2.6) 

is modular invariant. The easiest examples of modular invariants are M = I, and M = C, 
where C denotes charge conjugation. More generally, any symmetry of the unextended 
Dynkin diagram defines a modular invariant; on the other hand, symmetries of the extended 
Dynkin diagram (i.e. simple currents) may or may not yield a modular invariant for a given 
fixed level k (for example k must be even for su(2)). The symmetries of the unextended 
Dynkin diagram are often called conjugations. 

For a given modular invariant M, we call the exponents of M the diagonal terms M\\, 
or more precisely, the 'multi-set' Em in which M\\ copies of primary field A appears, for 
each A G P+. 

We shall be interested in constructing boundary states for these conformal field the- 
ories. A boundary state ||a)) is a coherent state in the full conformal field theory. In the 
simplest case where the boundary preserves the full affine algebra (that is generated by 
J^J it is characterised by the 'gluing condition' 

(Jm + J- m ) \\a)) = for all b and m G Z . (2.7) 



Since the modes that appear in ( [2.71) map each H\ (g) 7i* into itself, we can solve the 
gluing condition separately for each summand in (|2.5|). We can find a non-trivial solution 



provided that H\ is isomorphic to 7i M , i.e. for each A G Em- For each such A, the relevant 
coherent state |A)) is then unique up to normalisation, and is called the Ishibashi state. We 
shall always (partially) fix the normalisation of the Ishibashi states so that 

«A| q^o+Lo-^) ^ = S x , Xx (t) , q = e 2 ™ . (2.8) 
Every boundary state can then be written as a linear combination of Ishibashi states, 



V b 0fj, 



where we have introduced the factor involving Sq^ for future convenience. Given the above 
normalisation of the Ishibashi states, the boundary states are determined in terms of the 
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matrix ij) a n- In general we expect to be able to find as many boundary states as there are 
Ishibashi states, i.e. that ip is a square matrix. Furthermore, the completeness argument 
of suggests that i/j should be unitary. 

Not every linear combination (|2.9| ) defines an actual boundary state. The allowed 
boundary states have to satisfy a number of consistency conditions, the most important of 
which is the so-called Cardy condition jl| which can be understood as follow. The 'overlap' 
between two boundary states can be calculated from (|2.9|) as 



mq t(Lo + Lo-^) lh))= J2 %^X,(r). (2.10) 



Upon the modular transformation r i— > — 1/r, this amplitude must then be expressible in 
terms of a non-negative integer combination of characters, i.e. 

E %^X,(-Vr)= £ £ %^^XA(r)=: £ A/^(t), (2.11) 
where A/^ b , defined by 

JV& = £ ^f"^* . (2.12) 

fJ-EE M 

must be a matrix of non-negative integers. The collection of these matrices is usually called 
a NIM-rep. More abstractly, we define a NIM-rep A/" to be an assignment of a matrix A/a, 
with non-negative integer entries, to each A G P+ such that A/" forms a representation of 
the fusion ring 

A/aA/^ = £ iV^A/l, (2.13) 
for all primaries A, G P+ . Furthermore we require that 

(2.14) 

A/ca =A/a AgP+. 



' Actually, ( |2.8| ) only determines the normalisation of the Ishibashi state up to an arbitrary 
phase. Furthermore, if some \x appears with non-trivial multiplicity n in £m, we can redefine the 
Ishibashi states by a transformation in U(n). Unless these choices are specified further, there is 
therefore an ambiguity in the definition of if). This accounts for the difference in our formulae for 
ijj from those in flU 
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If the NIM-rep arises as in the above construction of boundary states with ifi being unitary, 
it is easy to see that ( |2.13|) and Q2.14j) follow from (|2.12|) . The easiest example of a NIM- 



rep are the fusion matrices themselves since the assignment A i— > N\ clearly satisfies both 
(fTTJD and ( CT) . 



The matrices Af\ of any NIM-rep can always be simultaneously diagonalised, by a 
unitary matrix ijj, in such a way that 

The sum will be over some 'multi-subset' £ = £{M) of P+, i.e. each element of P+ will 
come with a multiplicity, possibly zero. This multi-set £ is called the exponents of the 
NIM-rep M. The eigenvalues and simultaneous eigenvectors {i.e. the columns of ip) of 
the matrices M\ are parametrised naturally by the exponents \i G £{M). The entries of 
these eigenvectors, i.e. the rows of ip, are parametrised by x G £>; the same labels also run 
through the rows and columns of the matrices J\f\. 

As we shall see later on, some of the NIM-reps that we shall discuss in this paper have 
an analogue of a simple current symmetry (compare ( |2.4j) ). By this we mean permutations 
J B of B and J E of £, such that 

L^=e 2ligE(o) i 



an 



(2.16) 

Af b Xa ? =► Q E (X) + Q E (a) = Q E (b) (mod 1) , 



\fJ a b _ K rb 



for appropriate phases Q B and Q E . 

Any modular invariant (boundary) conformal field theory has both a modular invariant 
M and a NIM-rep N '. We say that the set of boundary states is complete, if M and M 
are compatible in the sense that their exponents agree (including multiplicities): £m = 
£{M). If this is the case we can then also identify the matrices ifi = this is to say, the 
diagonalising matrix ip of the abstract NIM-rep has then an interpretation in terms of the 
boundary states of the conformal field theory. 

One example where this identification has been understood is the diagonal modular 
invariant theory M = I for which £m = P+- The boundary states of this theory have 
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been constructed in |lj]; they are labelled by elements in P + and denned by ip = S. The 
corresponding NIM-rep is then the fusion matrix, as follows from (|2.2|) , which again has 
£{N)=P + . 

However, in general it is not known how to construct the NIM-rep that corresponds to 
the other modular invariants. In this paper we want to solve this problem, in the context 
of the affine algebras, for those modular invariants that correspond to a symmetry of the 
unextended Dynkin diagram. We will also get a natural interpretation for the 'boundary 
state labels' x G B. 



3. The abstract construction 

As we mentioned in the previous subsection, the boundary states for the diagonal modular 
invariant are known, and they are given by 

« A »= E -7r=lA*»- C 3 - 1 ) 

Suppose uj is any outer automorphism of q that arises from an automorphism of the 
unextended Dynkin diagram. We want to find the boundary states of the modular invariant 
based on u. The exponents S u of the corresponding NIM-rep are precisely the subset of 
P+ consisting of those representations that are invariant under the action of u. The 
problem of finding the boundary states for the modular invariant based on uj is therefore 
essentially equivalent to finding the w-twisted boundary states for the original diagonal 
modular invariant theory. A boundary state is w-twisted if it satisfies the twisted gluing 
condition 

{J b m + u (j*L m )) ||a)r = for all b and m G Z . (3.2) 

Twisted boundary states have been discussed before in [p3| . We shall propose a simple 
description for these boundary states (and their corresponding NIM-rep) below. 

Every twisted boundary state can be written as a linear superposition of the w-twisted 
Ishibashi states \fx)) w , that can be defined for each fx G £ w , i.e. for each \i G P+ that is 
invariant under the induced action of u, u* (fx) = fx. Again, the Ishibashi states are uniquely 
determined by this condition up to normalisation, which we (partially) fix by demanding 
that 

"«A| q hto+Lo-& \fx)r = 5 X , X x(r) . (3.3) 
10 



Let us expand the w-twisted boundary states as 

N>r= E -7r=l/*» w . (3- 4 ) 

where x is, at this stage, an abstract label for the different twisted boundary states. The 
NIM-rep that is associated to the cu-twisted boundary states is then simply given by 

A/L = E , (3.5) 



as follows from the same calculation that led to ( [2.12| ). 



Now we come to describe the main idea of our construction. The advantage of considering 
the problem in the context of the diagonal modular invariant theory (rather than as the 
problem of finding the boundary states corresponding to a different theory) is that we do 
not only have the Cardy condition that requires that the boundary states ( |3.4[ ) among 
themselves lead to non- negative integers as in ( |3.5|) . If the twisted boundary states (|3.4| ) 
are consistent, they must also satisfy an appropriate Cardy condition involving the overlap 
between one of the twisted boundary states, and one of the original Cardy states (|3.1|). As 
we shall explain momentarily, this constraint turns out to be very powerful, and will allow 
us to make a very suggestive proposal for 

In order to analyse this constraint we have to determine the overlap between the 
boundary state ||A)) and the boundary state Given our expansion in terms of the 

relevant Ishibashi states this can be easily done, and we find 

«M q^+Lo-^) || A)) = £ %^ x H( r) , (3.6) 
where X^\ T ) is the twining character in the representation //, i.e. 

Here t u is the induced action of ui on 'H VL . 

We can now use the fact that the above twining character agrees precisely with the 
ordinary character of the so-called orbit Lie algebra q J24j. (Explicit descriptions for these 



Lie algebras are given in table (2.24) of ||24| ; see also the translation of their notation to 



the notation of Kac J25| on page 12.) 
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In order to relate the 'closed string result' (3J5) to the 'open string picture' we have 



to apply again a modular transformation, as before in ( |2.11| ). In the present context 



this involves the characters of g that do not transform into one another under the S- 
modular transformation. Instead, the characters of g transform into linear combinations 
of characters of the twisted Lie algebra g w . Let us label the representations of g w by I. 
Then the relevant formula is 



xS" ) (-Vr) = E%xr(r), (3.* 



where S is the corresponding (unitary) S'-matrix that is given in Thms.13.8, 13.9 of |p5|| 
(except that for reasons of convenience we have chosen ours to be the transpose of Kac's). 
Every (ordinary) representation of the twisted algebra g w is, by definition, the same as 
a w-twisted representation of the original algebra g. Thus I labels equally the w-twisted 
representations of g (see for example |26| for an introduction into these matters). 



Putting all of this together, we can therefore write the overlap ( |3.6| ) in the 'open string 
picture' as 

2^ cT ■ ( 3 - 9 ) 

The relevant 'Cardy' condition is then that 

defines non-negative integers. 

Next we want to make a suggestive proposal for how to solve the two 'Cardy' conditions 
( pT5| ) and ( |3.10| ). The NIM-rep we are looking for has dimension |£ w |; because S is an 
invertible matrix, this equals the number of w-twisted representations of the conformal 
field theory. Now we can always construct a NIM-rep of this dimension by taking B to 
label the different u;-twisted (irreducible highest weight) representations, and by taking 
Af% x to be the twisted fusion rules that describe the fusion of the untwisted representation 
A with the twisted representation x, leading to twisted representations described by y. 



(The fusion of twisted representations is for example discussed in [p7j.) These twisted 
fusion rules automatically define a NIM-rep since fusion is associative. 

This suggests therefore that the label x that describes the twisted boundary states 
can be identified with the label for the w-twisted representations. Furthermore, the two 
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sets of integers ( |3.5|) and ( |3.10| ) should both simply be the twisted fusion rules. * Since 
these fusion rules are symmetric we therefore have, writing I = x, 

Comparing (|3.5|) and ( |3.10| ) this then implies that 

ifrxn = S xll . (3.12) 

Thus we propose that the boundary states of the NIM-rep associated to uj are labelled by 
the u-twisted representations, and the NIM-rep itself is precisely described by the twisted 
fusion rules. The matrix diagonalising the matrices Af\ is the modular matrix S. 

Returning to the original problem of finding the boundary states for the conjugation 
modular invariant, we can now write down a complete set of boundary states, all of which 
satisfy the gluing condition ( |2.7|) . As above, these boundary states are labelled by the 
a;-twisted representations of g, and they are given as 

where \fj)) are the standard Ishibashi states that exist for \i e S w . The Cardy condition for 
these boundary states is now precisely the condition that ( |3.5| ) are non-negative integers. 

Finally, we should mention that the above arguments imply that we have a generalised 
Verlinde formula 



= E — c*" S " - ( 3 - 14 ) 
~7 Son 



that describes the twisted fusion rules in terms of the S'-matrices and S Xil . 

3.1. Generalised fusion algebras 

It is natural to think of our construction in terms of some generalised fusion algebra 
that combines untwisted and twisted representations. Let us concentrate on the case where 



£ After this paper was completed we learned that the possibility of this interpretation in terms 



of twisted fusion rules was made in passing in [1C]. 
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the twist is of order two in the following. The fusion matrix corresponding to an untwisted 
representation A then takes the form 

Nr= ( N t A)- (3 ' 15) 

while the fusion matrix for a twisted representation x is 

n ™=(k % I'")' <3 ' 16) 

In both cases we have written the matrices in blocks corresponding to untwisted (A, v G 
P+) and twisted (x,y G <6) representations. These fusion rules define a consistent fusion 
algebra in the sense of || (although, unlike the more familiar fusion rings of conformal 
field theory, it is not self-dual). 

The NIM-rep that is defined by combining the untwisted and the twisted boundary 
states is precisely this generalised fusion algebra. In particular, this implies that the NIM- 
reps corresponding to M = I and the one corresponding to the modular invariant associ- 
ated to the conjugation uj combine into some larger algebraic structure. This assumption 
played an essential role in our 'derivation' above. 

This is quite analogous to the manner in which group theorists treat projective repre- 
sentations of a finite group G: they are interpreted as true representations of an extension 
H of the Schur multiplier M(G) of G by G. It would be interesting to try to push this 
analogy further and see, for example, what the Schur multiplier corresponds to in our case. 

Since the algebra defined by the matrices (|3.15| ) and ( |3.16|) defines a fusion algebra in 
the sense of [ffl, it possesses a (non-symmetric) S'-matrix that diagonalises the fusion rules, 
and recovers the full fusion coefficients iV fu11 via Verlinde's formula. In the present context 
this (\P+\ + \£u\) x (l-P+l + \£u>\) S'-matrix is simply given by 

where the first 2\E UJ \ columns are parametrised by \x G E u and either choice of sign (so each 
entry '^^m' appears twice), and where the final \P+ \ — \£ u \ columns are parametrised by 
the j G P+ \£oj. The first rows are parametrised by i G P+, and the last \£ w \ rows are 
parametrised by the w-twisted representations x of q. If we introduce the column vectors 
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then we can rephrase this as follows: for each \i and choice of sign, the column vector 
S Uv,±) is an ei § envector of N t n ( A G p +) with eigenvalue §^ = <S^ |±) /S^ i±) , and 
an eigenvector of A^ u11 (x G B) with eigenvalue ±-§^- = ^"j^ ±)/^o U (Jt ±)- Similarly, the 
column vector S^j is an eigenvector of iV| u11 (A G P + ) with eigenvalue = S^ l /S^ u , 
and an eigenvector of AT* u11 (x G £>) with eigenvalue = S^/Sqj 1 . 

This generalised fusion algebra also appears in || under the name 'classifying algebra'. 
Their equations (4.5) and (4.7) correspond here to ( |3.17| ) and (|3.15| ), (|3.16| ), respectively. 



One consequence of (|3.17|) and Perron-Frobenius theory is that, if the coefficients A/T 
are all to be non-negative, the entries S x o must be of constant sign for all x G B. We will 
see in the next section that they are indeed always positive. This can be regarded as a 
non-trivial test of non-negativity — the only property of our NIM-reps which we cannot 
prove in general. The property of S x q to be positive has also a physical interpretation since 
ifixo = S x o is proportional to the (positive) boundary entropy of the boundary condition 
described by x [ P8| . 

3.2. Generalised NIM-reps 

Defect lines (that is, non-local operators attached to non-contractible loops) can be 
inserted into the partition functions of the torus and cylinder, resulting in two other sets of 
NIM-rep-like data associated to a conformal field theory |]] . These also appear naturally in 
the subfactor context (see for example Thm. 4.16 in ||16|| ). In particular, to each pair A, A' 
of weights in P + , we should have a matrix V\\> = (^AVi)ijeB °^ non-negative integers, 
giving a representation of the 'double-fusion algebra' 

v X yV^= J2 Wy^> ( 3 - 19 ) 

and obeying both V^ cx ,.i = M\y and V£ x , = Vc\,c\'- Here M is the torus partition 
function of the theory, C denotes charge conjugation, and V t is the transpose of V . The 
label 1 is a special element of B. These conditions imply that the matrices Vw can be 
simultaneously diagonalised by some matrix -0, and that the eigenvalues of V\\> will be 

5* S ! ! 

So s ^ ' ^ S mu hiplicity is required to be (M^^/) 2 . The other NIM-rep-like data concerns 
a representation by non-negative integer matrices Mi = (A/^)x,ye6, i £ B, of a (not 
necessarily commutative) fusion-like ring built out of the matrix ip — see section 4.2 of 
for details. 
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In the case where M is a permutation matrix M\n = S^^x (this is, in particular, the 
case when M is a conjugation modular invariant, as considered in this section, or the simple 
current modular invariant considered in section 6.1 below), it is easy to find these matrices 
V\\> and Hi- in particular, we can naturally identify B with P + , and put Vw = NxN-^x' 
and Af\ = Af\. Thus these other NIM-rep-like data can also be constructed for the theories 
considered in this paper. 



4. Explicit formulae 

In this section we shall collect explicit formulae describing the NIM-reps for the various 
classes of algebras. The NIM-rep is uniquely determined once we know the exponents £{M), 
the eigenvalues for all weights A G P+ and exponents \i G £(AT), and the diagonalising 
matrix ip x ^- As we have explained above, the matrix ijj can be identified with the modular 



S'-matrix of the twisted algebra g 10 . In principle these 5"- matrices are known |p5[ , but 
as for the case of the more familiar untwisted S'-matrices, the usual formulae are rather 
complicated (involving sums over Weyl groups). In the following we shall therefore give 
simple expressions for both S'-matrices, following the ideas of ||29|1 ■ We shall also explain 
how the NIM-rep coefficients can be expressed in terms of the usual fusion coefficients; 
this will make manifest that they are integers, as they must be. We conjecture that they 
are also non-negative, and have verified this on the computer for several algebras and 
levels. The arguments behind these calculations are given in the following section, and 
some fusion graphs are exhibited in the appendix. In the following we shall analyse the 
different algebras case by case; we shall always use the same numbering of nodes as in 



4-1. The A-series 

A level k weight A G P + for the A-series, i.e. su(n) = A n _i, looks like A = J^™^ 1 AjAj 
where ^ Aj = k. Charge conjugation 10 = C for su(n) takes the weight A to the weight 
CX = AoAo + J27=i ^n-i-h-i- Charge conjugation is non-trivial for su(n), except for su(2). 

The exponents S here are the C-invariant weights, i.e. those /i with Dynkin labels 
[i n -i = li% for 1 < i < n — 1. All of these appear with multiplicity 1. The eigenvalues, 
involving ratios of S'-matrix elements, can be effectively computed as follows. 
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To the weight A = J2i K-h-i associate orthogonal coordinates X + [£] = n — j + Yl7=e ^* 



for l = \, 



n, 



and write t + (X) 



n (n — 1) 
2 



s 



x^ = s exp 



2ni 



(k + n)n 



det I exp 



-2yri 



Then 

A + [i]^ + [j] 
k + n 



(4.1) 



l<i,j'<n 



where s is some (for our purposes irrelevant) constant that is given in e.g. H|. Here 'det' 
denotes the determinant of the n x n matrix whose (i, j)-th entry is provided. 

In order to describe the formula for ip, we need to distinguish two cases depending on 
the parity of n. 

The case of su(2n + 1) 

As we know, the exponents, which parametrise the eigenvectors, and label the columns 
of if), are the C-invariant weights of q = su(2n + 1) level k. Alternatively, we can think 

(2) 

of them as the level k weights of the orbit Lie algebra q = A 2r f. In either case, they 
can be equated with the (n+1) -tuples \x = (ho; Hii • ■ ■ , A* n ), where \ii G Z>o and we have 
k = no + 2 hi + ■ ■ ■ + 2Hn-i + 2Hn- To each exponent /U E £ we associate the coordinates 



H[i\ 



n 



+ 1 — i + YTj=i H for i = 1, . . . 



n. 



The boundary states, which label the components of each eigenvector, and both the 
rows and columns of each NIM-rep matrix Af\ = (A/^ a ), as well as the rows of if), are the 
level k C-twisted weights of q = su(2n + 1), or alternatively the level k weights of the 

(2) 

twisted Lie algebra q w = A\^. They can be equated with all (n + l)-tuples (ao; ai, . . . , a n ) 
where k = ao + 2ai + 2a2 + • • • + 2a n and a, G Z>q. To each such a we associate the 



coordinates a\i] 



n 



+ l — i + J2™=i a j for i = 1 



n. 



Then the ip-m&trix, that is to say the modular matrix S for A 2r f, is 



(2) • t 



a jU 



n (n — 1) 



(k + 2n+ 1)2 



det 



sin 



2na[i\ h[j\ 
k + 2n + l 



l<i,j'<n 



(4.2) 
0), and 



We will show in section 5 that the column of if) corresponding to h° = (k;Q, 
the row of if) corresponding to a = (k; 0, . . . , 0) are both strictly positive. 

Incidentally, this ^-matrix can be regarded as a symmetric submatrix of the S'-matrix 
for B n level k + 2, provided we identify a and h with the appropriate weights a', h' m 



* We shall always use the convention that =0. 

t This determinant formula for S for the case of first appeared in [24|. They also noted 



that when k is odd, S equals the 5-matrix for C n level ^2^"- 
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P+ +2 (B n ). More specifically, we identify p \— > p' = (po + /-*i + 1; pi, ■ ■ ■ , p n -i, + 1) an d 
a i— > a' = (ao + ai + 1; 01, . . . , a n _i, 2a n + l).* With these identifications, the ^-matrix 
becomes symmetric, and the relation between ip and the B n level k + 2 S'-matrix 5" is 
simply W =2S' a>fM ,. 

Using this identification, we can express the NIM-rep coefficients A/^ a , corresponding 
via (|3.5|) to (|4.2p ) in terms of the ordinary fusion coefficients N' of B n level fc + 2, and 
branching rules for B n C su(2n + 1), 

*L = E 6 " (KJ' - N JY,a' b ') ■ (4-3) 
7' 

Here, J' is the simple current of B n , which acts on weights v' by Ju' = [y' x ; v' , u' 2l . . . , 
The coefficients by describe the branching rules ©y 6y (7') = (A) for the embedding of 



B n C su(2n + 1); some branching rules for B n C su(2?t, + 1) are given for instance in pi 



The sum in ( f4.3| ) is over all dominant weights 7' of B n ; we extend the definition of the 



fusion coefficient NL a , b to arbitrary dominant weight 7' in the obvious way (compare the 



7 ! a 



Kac- Walton fusion formula): the fusion coefficient is det(io) , a , b if there is a unique 
element w of the affine Weyl group of B n for which 10.7' := 10(7' + p) — p lies in P+ +2 (B n ), 
and it is zero if there are more than one such w's (i.e. if 7' lies on the boundary of the 
Weyl alcove). Equation ( |4.3| ) makes it manifest that our coefficients N\ a are integers. 



When the level fc is odd, a simplification occurs in that we can identify the boundary 
state a with the entire C n level weight a = ( a °~ 1 ; ai, . . . , a n ), and similarly for p. 
With this, ip becomes the C n level S'-matrix, and the NIM-rep becomes 



N\ a = N ld h , (4.4) 



Z^7 



where ©6~ (7) = (A) are the branching rules for C n C su(2n + 1). As for (|4.3|) , the fusion 
coefficients N^a b are defined for 7 not necessarily in P + , in which case they can be negative. 
Some branching rules for C n C su(2n + 1) are given in |^l] . 

This simplification reveals a hidden symmetry in the NIM-reps at odd level fc: the C n 
simple current J =: J B acts on a E B = P + 2 (C n ) by J B a = (2a n + 1; a n _i, . . . , ai, a ° 2 " 1 ) 
with phase Q s (^) = § ELi *A*i, and A/"£X = A^ a (see flOED ). 

Some low-rank clarifications are needed in both Ql.3| ) and ( |4.4| ). By I?i level fc + 2 
we mean affine su(2) at level 2(fc + 2). Also, su(2)=so(3) is embedded in su(3) in two 



P+(B n ) consists of the weights p! = (p' ; p[ , . . . , p' n ) for which p' + p[+2 5^ n =2 Pi+Pn = k. 
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different ways; one has branching rule (2) = (1,0) and the other has (0) © (1) = (1,0). 
By Bi C su(3) we mean the former embedding, and by C\ C su(3) we mean the latter. 
Similarly sp(4)=so(5) is embedded in su(5) in two different ways: one with branching rule 
(1,0) = (1,0,0,0) (using the so(5) Dynkin labelling) and the other with (1,0) © (0,0) = 
(1,0,0,0) (using the sp(4) Dynkin labelling). By B2 C su(5) we mean the former and by 
C2 C su(5) we mean the latter. 

The case of su(2n) 

Next we determine the formula for ijj for g = su(2n) at level k. The twisted algebra g^ 

(2) (2) 
here is g w = A 2r l_ 1 (for n > 2), or with nodes 1 and 2 interchanged (for n = 2). 

The rows of ifi are now labelled by (n + l)-tuples a with k = cio + ai + 2a^ + ■ ■ ■ + 2a n . 

To each such a we associate the coordinates a[i] = n + 1 — i + Yll=i a j for z = 1, . . . , 71. 

(21 

The orbit Lie algebra is g = D^!^, each exponent [i G £ can be regarded as an (n+1)- 
tuple /i, where k = /j,q + 2\i\ + ■ ■ ■ + 2\i n _\ + \i n . To it we associate the coordinates 
= 2n + 1 — 2i + 2 5^<=j + £*n- Then the V'-matrix (i.e. S for A 2r l-x) * s 

2 

^ afl = (-l) 1 ^ (k + 2n) _ * 2 n "5 det 

Although ift is non-symmetric, both the column corresponding to fjP = (k; 0, ... , 0) and 
the row corresponding to a = (k; 0, . . . , 0) are strictly positive. 

We can interpret the resulting NIM-rep coefficients J\f\ a using fusions N' and branch- 
ing coefficients b^, for the semisimple algebra su(2)ffi- • -®su(2) (n times) at level k + 2n — 2 

d 7' 7T 

This formula will be explained in more detail in section 5. We can replace the sum over 
7T with the determinant of an n x n matrix made up of various A% fusions. Even so, 
equation (|4.6|) is not very practical, but it does make it manifest that our coefficients Af^ a 
are integers. 

(2) 

The symmetry of the Dynkin diagram of g = D^+x gives rise to a simple current 
symmetry J E of £ , while the symmetry of the Dynkin diagram of g w = A^—i yields a 
simple current symmetry J B of B — see ( |2.16| ). These will be described in section 5.3. 
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. ( ir a[i] n\j]\ 
sin - — — 
V k + 2n J 



4-2. The D-series 



Next we consider the algebra sb(2n) = D n at level k. A weight A G P+ here satisfies 
k = X Q + Ai + A n _i + A n + 2 Xir=2 2 A;- 

As before, we need to compute effectively certain ratios of S'-matrix elements. In order 
to do so we associate to the weight A = (Ao; Ai, A2, • • • , A n ) the orthogonal coordinates 
X+\i\ = n - £ + Y?l} A, + A "~ A - 1 for I = 1, . . . , n. Then 



S 



det 



cos ( 2n 



A + [i]/x+[j] 
k + 2n - 2 



+ (-i) n det 



1<z,j <n 



sin ( 2n 



A + [»]/i + [y] 
k + 2n - 2 



l<i,j'<n 

(4.7) 



where s is some irrelevant constant (which is again given in e.g. ||29|| ). 



The case of so(2n) with chirality flip 

The chirality flip uj = P interchanges the nth and (n — l)st Dynkin labels, A n <-> A n -i- 

The chirality flip P agrees with charge conjugation C = S 2 when n is odd, but C is trivial 

(2) 

when n is even. The relevant orbit Lie algebra here is q = ^2n-3 - ^he exponents are 
the P-invariant D n weights, i.e. the weights of the form (//q; jUi, A*2 5 • • • 5 Hn-x, l^n-x) where 
= ^0 + ^1 + 2^2 + - • - + 2fx n -i. Given any P-invariant Dn^ weight \i, define its coordinates 
by 



n 



i + EL/ ^j, for i = 1, . . . , n - 1. 



(2) • , 
The relevant twisted algebra here is = D n . The boundary states are all n-tuples 

(ao; ai, . . . , a n -i) °f non-negative integers such that k = ao + 2ai + 2a2 + - ■ •+2a n —2 + a n -i- 
Define coordinates a\i\ = 2n — 1 — 2i + 2T]™ = f 



■^-matrix (i.e. the matrix S of D n ) is given by 
1 



1 — 2i + 2 Y^j=i a j + a n-i fori = 1, . . . , n — 1. Then the 

,(2)> 



'0 



fc + 2n - 2 



rt-l 
2 



2 n-2 ^2 i« +"+ 2 det 



sin 



k + 2n-2 



l<i,j'<n— 1 



Again %j) is non-symmetric, but both its fjP = (k; 0, . . . , 0) column and a = (A;; 0, ... , 0) 
row are strictly positive. 

The NIM-rep coefficients can be interpreted in terms of fusions N' of P n _i at level 
k + 1, and branching coefficients 6^, of P n _i C so(2n), 



(4.9) 



This will be explained in section 5.2; as before it demonstrates that our coefficients M\ a 
are integers. 
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(2) 

The symmetry of the Dynkin diagram of g = ^4 2 n_i gi ves rise to a simple current 
symmetry J E of £ , while the symmetry of the Dynkin diagram of g u = Dn yields a 



simple current symmetry J of B — see ( |2.16| ). These will be described in section 5.3. 



The case of so(8) with triality 

There are 4 other non-trivial conjugations for so(8), but their NIM-reps are determined (see 
the end of this subsection) once we know the NIM-rep corresponding to 'triality' uj = T. 

The order-3 triality maps the Dynkin labels fx = ((Xq', /-ti, /U 2 , /U3, ^4) to (^0; A*4? ^2, f^i, 1^3)- 

(3) ^ 
The relevant twisted algebra here is g = g u = D\ . The exponents are the D4 weights of 

the form (/i ; fix, /U 2 , A*i, A*i) where k = liq + 3/ii + 2/i 2 . Given any such exponent //, define 

coordinates fi[l] = 3/Ui + ,u 2 + 4, «[2] = ,u 2 + 1- 

The boundary states are labelled by all triples (ao; 01,02) of non-negative integers 

such that k = ao + 2a± + 3a 2 . Define coordinates a[l] = a\ + a 2 + 2, a[2] = 02 + 1. Put 

k = + 6. Then the ^-matrix (i.e. the matrix S for -D^ 3 " 1 ) is 

=— {ciram! + 2mn + 2nm + nn') + c(2mm! + ran + nm' — nn') 

K V 

+ c(— mm' + mn' + nm' + 2nn) — c{2mm + ran + nm' + 2nn) (4-10) 
— c(mm' + 2mn' — nm' + nn') — c{mm! — mn + 2nm + nn')j , 

where c(x) = cos(^^), and we take m = a[l], n = a [2], m' = u[l], n' = fi[2j. Again tp 
is non-symmetric, but both its fx = (k; 0, 0, 0, 0) column and a = (k; 0, 0, 0, 0) row are 
strictly positive. 

We can find an interpretation of the NIM-rep coefficients M in terms of fusions N' 
and branching coefficients iA of A 2 at level k + 3 

2 

•ML = EE 6 " (^W,</ - ivs, VjC , a , 6 ') . (4.11) 

i=0 7' 



This is further explained in section 5. The formula (|4.11|) implies, in particular, that our 
NIM-rep is manifestly integral. 

We can now give the NIM-reps for the 3 remaining conjugations of so(8). The NIM- 
rep for the conjugation T _1 : Ai — ► A4 — ► A3 — > Ai equals that for triality T, given above. 
The NIM-rep for the conjugation T _1 P : Ai *-> A3 is given by A 1— > Af<f x , where Af p is the 
NIM-rep for chirality flip, while that for TP : Ai <->• A4 is A 1— ► Af^-iy 
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4-3. The algebra Eq 



The final algebra with a non-trivial conjugation is Eq. Its level k weights A satisfy k = 
Ao + Ai + 2X2 + 3A3 + 2A4 + A5 + 2A6- The order-2 charge conjugation u = C interchanges 
the 1st and 5th, and 2nd and 4th, Dynkin labels, and fixes the Oth, 3rd, and 6th. The 
relevant twisted algebra here is g = g w = Eq . 

The ratios of S'-matrix elements can be effectively calculated using the formula 
of ||2g|| . The exponents are the C-invariant Eq weights, i.e. the \i of the form 
(fio] A*i, A*2, ^3, A*2, Mi, Me) where k = (Iq + 2\i\ + 4(i2 + 3(i3 + 2\iq. Define coordinates 
by (i[l] = 2/Lti + 3^2 + 2fi 3 + fi 6 + 8, (i[2\ = [i 2 + (13 + (i 6 + 3, //[3] = (i 2 + (13 + 2, and 
fl[4] = (Ji 2 + 1. 

The boundary states are all quintuples (ao; ai, 02, 03, 04) of non-negative integers such 
that k = ao + 2a\ + 3a 2 + Aa 3 + 2a<±. Define coordinates a[l] = a\ + |a2 + 2a3 + 04 + 



a [2] = \a 2 + a 3 + a 4 + |, a [3] 
the ^-matrix (i.e. the matrix 5 1 for E ( 

'2na[{\ 



\a,2 + 03 + f , and a [4] 
is 



±a 2 + h. Put k = k + 12. Then 



(2), 
6 , 



det 



sin 



K 



+ det 



l<i,i<4 



sin 



2tt (c.a)[i]/z[)'] 



l<i,j<4 



+ det 



sm 



2n(c t .a)[i\n\j) 



K 



(4.12) 



J l<i,j<4 j 

Here, (c.a) [i] denotes the ith coordinate of the matrix product of c with the column vector 
with entries a[j], and c* means the transpose (=inverse) of c, where c is the orthogonal 
matrix 

1 

c = - 



1 - X \ 

-1 1 

1 1 

-1 -1/ 



(4.13) 

0) row is strictly positive, and the same 



1 

1 1 
1 -1 
\l -1 

Again, ift is non-symmetric but the a = (k; 0, 
appears to hold for the (i° = (k; 0, ... , 0) column. The positivity of the a row can be 
shown in the usual manner, but the arguments to prove the positivity of the n column 
are much messier; we have therefore only verified this claim on a computer for levels up to 
k = 100. 

We can interpret our Eq NIM-reps using so (8) fusion coefficients N' and branching 



coefficients by 



J> 7' 7T 



A a 



b' 



(4.14) 



This sum (|4.14| ) is manifestly an integer. 
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5. Calculations and further elaborations 



In this section we want to provide details of the derivation of the various formulae 
we gave in the previous section. First we shall explain how we obtained the formulae for 
S = ip. As an example we shall discuss the calculation for the algebra g = su(2n) explicitly. 

5.1. Derivation of the ^-formula for g = su(2?i) 

From the discussion of section 3, we need to identify the quantities of Thm. 13.9 in 
PS[] . The finite Weyl group of both g — D^-y and g^ = A^-i ^ s J us t the Weyl group 



of sp(2n), so the ifi formula will involve a determinant that will look like the one in the 
S'-matrix for sp(2n). Furthermore, we need to be able to evaluate inner products (a\fx), or 
equivalently we need the values of (Aj|A^) for the fundamental weights Aj,A^ of g w and 
0, respectively. To do this we use the equation (see page 222 of fl2"5fl ) that (aJ|Aj) = dij 
for the coroots ctj of g w , where the values of (ct/|aj) are determined from the Cartan 
matrix of ^n-i ( see (6-2.1) of pq] ). We obtain the realisation Aj = e\ + ■ ■ ■ + ei for all 
i = l,...,n, where is an orthonormal basis of IR n . We now find A^ = A^ for i < n, 
and A^ = |A n , by page 266 of f25[l . This implies that we should define the orthogonal 



coordinates fi[i] and a[i], as we did in subsection 4.1, and rescale the inner products by 2. 

The overall multiplicative factor is computed as follows. The number of positive roots 
of the horizontal algebra C n of is n 2 . The lattice M' is the lattice spanned by the 
coroots a^, which is the square lattice Z n . The lattice M is the lattice spanned by the 
roots «i, which is the D n root lattice (see page 93 of H). The index \M'/M\ is y/2, and 



the index \M*/(k + h v )M\ is 4 (k + 2n) n . We want to replace e lz — e lz with — 2isin(;s), 
and so we get an additional factor of 2 n i~ n . Collecting all these terms, we recover (^ 

5.2. Relation of NIM-reps to ordinary fusion rules 

Now we turn to the derivations of the fusion formulae for our NIM-reps. Consider 
first the easiest case: su(2n + 1) when k is odd. Here, we noticed that ifj a ^ = S'a^, where 
the tilde denotes quantities for C n level ^-k The map u i— > u described in section 4.1 
is onto all of P + 2 (C n ). We know that an su(2n + 1) character, restricted to the Cartan 
subalgebra (CSA) of C n C su(2n + 1), equals the sum of C n characters, the sum given by 
the appropriate branching coefficients 6~ . Now, the su(2n + 1) ratio is the A-character 
of su(2n + 1) evaluated at the element — 27ri k _^2n+i m the CSA not only of su(2?t, + 1), 
but in fact of the embedded algebra C n , since Cfi = \i. This means that it equals the sum 
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&~ Styp./Sfifi of C n characters evaluated at that element of the CSA of C n . Thus we 
have succeeded in expressing all su(2n + 1) quantities in (|3.5| ) in terms of C n quantities, 
and we obtain ( |4.4| ) from [5(J and 



where here (as throughout this paper) we extend the definition of fusion coefficients and 
S'-entries to arbitrary dominant weights, as explained after ( |4.3| ). 

Next we want to verify ( ^.3| ). The only modification that arises relative to the previous 
case is that the set of all \x' is only a subset of all weig hts in P' + := P+ +2 (B n ), namely all 
those weights with fx' n odd and fi' Q > \i' x . In order to recover the B n fusions from ( |3.5[ ), we 
need to extend the sum from all these \i' associated to n G E, to all weights // G P' + . This 
is where the simple current J' comes in. Its phase Q(fi') is ^f- (recall ( |2.3|) ), and so the 
alternating sum in ( |4.3|) projects away \j! n even, as well i , and absorbs an extra 

factor of 2 coming from our formula ijj = 25". The weights with fi' < \j! x are recovered by 
taking J'-orbits of //, and this absorbs the remaining factor of 2. We have also used here 
that Q{i) G Z and Q(a'),Q(V) G Z + ±. 

The case of su(2n) 

In this case, the ^-matrix is a submatrix of S (rescaled by \/2) for B n level k + 1, using the 
identifications \i h-> (^o + A*i + 1; A*i, • • • , V>n) £ P+ +1 (B> n ) and a h-> (a ; 2a n + 

1) G P + +1 (.B n ). But in order to interpret the resulting NIM-rep coefficients Af£ a using 
ordinary fusions N', we need to break the algebra down further. 

Write A"l for the direct sum su(2) © ■ ■ -®su(2) (n times), and put P| = Pj; +2n ~ 2 (Ai). 
Given an exponent \x G E, define the weight jj! = . . . , jj!^ n ') G P' + by //w = (/c + 2n — 

2 - /*[*]) e i+ +2n " 2 (^i)- Similarly, define a 1 G P| by a'« = (k + 2n-2- a[i];a[i]) G 
P + +2n_2 (j4i). The quantities by in ( |4.6| ) are the A™ C su(2n) branching coefficients 
©7' by (Y) = (A). The first sum in (fO|) is over all n-tuples d G 1^ with £\ dj even. Ai 
has a simple current J' which sends z/ to (z^; z/q), and by J /d 7' we mean the weight with 
(J' d 7 ')« = J'^(yM). The third sum in ( fOP is over all permutations 7r G 5 n , which act 
on a' by (7ra')W = a'^ % \ The quantity e(7r) is the sign of 7r. 



Expand out (|4.5|) into 



V / -TT 



7T Z 7T 
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The set of all // here consist precisely of all the weig hts in P' + with > ^'i (2) > • • • > 
and = ■ ■ ■ = ii'x (mod 2). Again, we need to extend the sum over these jj! to 
a sum over all P\. The non-trivial simple current of Ai has phase Q{y) = Thus the 
sum in (|4.5|) over the 2 n_1 n-tuples d will retain only those \j! with /t^W = • • • = y!-^ 71 ' 
(mod 2), and with the numbers fi'^ all pairwise distinct. It also absorbs the extra 
factor of 2 n_1 coming from (5^2). An additional sum over the S'n-orbits of // removes the 
fi'i^ > fJi'i ^ > • ■ ■ > fJ'i^ restriction. This also cancels an extra factor of n\ coming from 
our over-counting of fusion coefficients due to the sum over ^-orbits of both a' and b' - 
we use here the fact that the branching coefficients by and b*, are equal, for any tv. 

The case of so(2n) 

Write P' + = P+ +1 (B n -i). Under the identification \i \— > jj! = (/xq; /■*!, • • • , fJ-n-2, 2 / u n _i + 



1) G P' + and oh>o'= (ao + «i + 1; «i, • • • , a n _i) G P+, the ^-matrix (|4.8|) is seen to be a 
submatrix (up to a factor of a/2) of the S'-matrix of P n _i at level k + 1. 

Recall that the simple current J of P n _i acts on 7' by interchanging 7q <-> 7^, and has 
phase QiY) = The // consist of all weights in P' + with fx' n _ 1 odd (z.e. the spinors). 

The coefficients by in (|4.9| ) are the branching rules Q)y e p^ by(Y) = (A) for the embedding 
P n _i C D n . Some simple branching rules are (A[) + (0') = (Ai), (A-) + (A-_ x ) = (A;) for 
1 < % < n — 1, and (A^_ x ) = (A n _i) = (A n ); others can be found for instance in 

The case of so(8) 

Using the identifications /i 1— > (/io! A*2> 3/ii + 2) and a 1— > (ao + ai + a2 + 2;a2,ai), the 
•0-matrix can be seen to be a submatrix of the level k + 2 S'-matrix of G2 (rescaled by \/3). 
It is clear from this interpretation of i/j that the a row of ijj will be strictly positive. That 
the fjP column, corresponding to the G2 weight (fc; 0, 2), of ip is strictly positive, is a little 
more awkward: we need to show 



c(2m + 3n) + c(3m + n) + c(m — 2n) > c(3m + 2n) + c(2m — n) + c(m + 3n) , (5.3) 

for all m = a± + ai + 2, n = + 1, where here c(x) = cos(^^); to do this, reduce it first 
to the quadratic inequality 

c{5) 2 - (2c(a) 3 - c{a)) c(5) + (4c(a) 4 - 3c(a) 2 ) > , (5.4) 

where a = m ^ 21 = " 1+2 2 a2+3 and 5 = 22Z y B = It can be shown by standard arguments 

that this final inequality is indeed satisfied for all a G B. 
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To get an interpretation of the NIM-rep coefficients M in terms of ordinary fusions 
N' , we must break G2 down to A 2 . In particular, write P' + = P^ +3 (A2). Associate the 
exponent fx G P+(D±) to the weight fx' = (/xq) \i 2 , 3/xi+/X2+3) G PL, and the boundary state 
a to the weight a' = (ao+ai+a2+2; a 2 , ai+a2 + l) G P' + . In ( |4.11[ ), J' is the simple current 
of A2, which sends v' to (u' 2 ; v'q, C there is charge conjugation, which interchanges 
p[ <-> z^ 2 - The quantities fey are the A2 C -D4 branching coefficients ©y fey (7') = (A). 
There are two of these branchings; the one we need comes from A2 C G2 C S3 C D4. 
Many of these branching rules can be found in ||31|| ; some of the more useful ones are 



2(0, 0) © (1, 0) © (0, 1) = (1, 0, 0, 0) = (0, 0, 1, 0) = (0, 0, 0, 1) and (1, 1) © 3(1, 0) © 3(0, 1) © 
2(0,0) = (0,1,0,0). 
Then ( [4.10] ) becomes 



expressing if) in terms of A2 data. 

The fx' consist of all weights in P' + with fx[ < fx' 2 and fx[ = fx' 2 (mod 3). The first 
condition on fx' is removed by the sum of C"-orbits of fx' - - we use the fact that the 
branching coefficients b\ and fejX, , are equal. The second condition on fx' is removed by 



The case of Eq 

In this case the ^-matrix is a submatrix of the level k + 3 S'-matrix (rescaled by 2) of 
F4, using the identifications fx 1— > (fio; fx^, fis, 1fX2 + 1, 2^ii + 1) and a t— > (ao + as + 0,2 + 
ai + 3; a4, as, 0,2, ai). Now, any F4 level k quantities can be rewritten in terms of sb(8) 
level k + 3 ones, using the embedding so(8) C F4, and this will enable us to interpret these 
NIM-reps using so(8) fusion coefficients. 

In particular, write P' + = P^ +6 (D4), and associate the exponent fx G P+(Eq) to the 
weight fx' = (fxo; 2fi\ + 2fX2 + fis + 4, fx§, fis, 2fi2 + fxs + 2) G P' + and the boundary state 
a to a' = (ao + ai + a2 + a^ + 3; ai + 02 + as + 2, a$, as, a2 + as + 1) G P' + . The sum in 
( fO^D over J' is over the 4 simple currents of so(8): the identity; J' v v' = (z/^; u' , v' 2 , v'^, u' 3 ); 
J' s v' — (v' li ;i' 3 ,i> 2 ,i'' l ,v' ); and J' v o J' s . The sum there over 7r is over all 6 conjugations 
of so(8), corresponding to each of the 6 possible permutations it of the 1st, 3rd, and 
4th Dynkin labels, and e(7r) is the sign of that permutation. The numbers fey are the 
-D4 C Eq branching coefficients ©y fey (7') = (A). Some useful D4 C Eq branching rules 



are 3(0000) © (1000) © (0010) © (0001) = (100000) = (000010) and 2(0000) © 2(1000) © 
(0100) © 2(0010) © 2(0001) = (000001). 




(5.5) 



the sum over i in (|4.11| ) (recall that the phase for J' is Q(v') 



3 



2 
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To prove ( [4.14|) , note first that if; a ^ = 2 e(n) S' na , where the sum is over all 6 



permutations tv of the labels a'i,a' 3 ,a' 4 . The weights jj! consist of all weights in P' + with 
Hi > h's > h'a and Hi = H3 = ^4 ( m od 2). The phases for the simple currents J' v and 
are Q v (v') = Vz 2 " 4 and Q 8 iy') = ~ 2 ^ 3 ■ The rest is as before. 

5.3. Further remarks 

It remains to show that, for the case of the A- and D-series, the column of ijj corre- 
sponding to fx = (k] 0, . . . , 0), and the row of ip corresponding to a = (k; 0, ... , 0) are 
both strictly positive. The key observation for the proof of these statements is the fact 
that ip is a submatrix of some untwisted S'-matrix (see section 5.2). In some cases (e.g. 
for the h° column for the case of su(2n)) strict positivity then follows automatically from 
the fact that the vacuum row or column of that S'-matrix is strictly positive (see Remark 
13.8 in |25|]). However, in most other cases we need to know the sign of the A n row and 
column of the ^-matrix for B n . 

Consider for concreteness the n° — (k; 0, . . . , 0) column of ip, for su(2n + 1). We know 
from section 4.1 that ^ aM o = 2S' a , A , . The sign of S' a , A , equals the sign of S' A , a , / 'S' , a ,, 
which is a value of the A n -character of B n . This is easy to calculate, and we find 



Given the definition of a[i], we therefore conclude that ip a 'u° > 0. All other algebras are 
handled in the same way. 

Finally, let us discuss the simple current symmetries for the various NIM-reps. The 
simple current J E for su(2n), which corresponds to the Dynkin symmetry of Dn \ permutes 
H e £ by J e h = (fx n ; (J>n-u ■ ■ ■ > A*o, fJ>u ■ ■ ■ 5 Vn-i)- Ks phase Q E (a) = \ J2i ia i wil1 be a 
grading for the NIM-rep coefficients (see ( |2.16|) ), where Q E (X) = i 5^i=i 1 The simple 
current J B for su(2n), which corresponds to the Dynkin symmetry of A^_ lt permutes 
a 6 B by ao «-> ai, and has phase Q b (h) = Hn- It defines a symmetry of the NIM-rep (see 

The simple currents J E and J B for so(2n) are defined similarly, except with their 
formulae interchanged (for example J E acts on \i by interchanging fiQ f^lj an d we have 
Q E (a) = a n ). The NIM-rep grading tfLWj) arises with the choice Q E {\) = ^(A n _i + A n ). 
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6. Generalisations 



The above constructions generalise also to certain other classes of modular invariants, in 
particular those that come from simple currents ( p.3|) . Here we shall only give the details 
for two simple complete description will be given elsewhere. Together, these two 

cases cover all su(p) level k, where p > 3 is prime. The 'i/'-matrices for all modular invariants 
of su(2) were given in |§ (i/j for the 'tadpole' su(2) NIM-rep, which does not correspond to 
a modular invariant, is a special case of the next subsection). The simple current NIM-reps 
given below are presumably consistent with those in [|11| , but we are more explicit here. We 
should also like to emphasise that, unlike us, [JTJ do not have a rigorous proof that their 
conjectured NIM-reps are non- negative integers (see the discussion in the introduction). 

Consider any su(n) level k. It has a simple current J of order n, corresponding to 
the cyclic symmetry of the extended Dynkin diagram, which sends the highest weight 
A = (Ao; Ai, . . . , X n -i) to JA = (A n _i; Ao, . . . , A n _2). For any divisor d of n, define the 
matrix M[d] by 

n / d / d ' k'\ 

M[d] x ,, = J2^ + <W>«a, (6-1) 

j=l \ / 

where t(X) — *Aj, and k' = k + n or k' = k depending on whether or not both k and n 
are odd, respectively. In ( |6.1[ ) we write S y (x) = or 1 depending, respectively, on whether 
or not - G Z. Then M[d] is a modular invariant if and only if the product (n — 1) kd is 
even. For instance, M[n] = I is a modular invariant for any su(n) level k; for su(2), M[l] 
is a modular invariant if and only if k is even. 

6.1. A case without fixed points 

Let us consider su(n) level k, for the simple current J, when gcd(/c,n) = 1. We will 
explicitly give below the corresponding -^-matrix and prove it yields a NIM-rep (as always, 
the hardest thing to prove is that the entries are non- negative integers). For such 
su(n) level k, the matrix M[l] in (|6.1| ) will be a modular invariant if and only if n is odd. 
However, the matrices defined by the -0-matrix given below define a NIM-rep even if this 
is not the case. For example, if n = 2 and k is odd, they describe the 'tadpole graph' of 
0. So for n even, this NIM-rep generalises the tadpole. 

Let Pq be all level k weights \x = (//i, . . . , \i n -\) of su(n), with n dividing t(n) = J2i Wi- 
Then Po will label both the rows and columns of ip. (We see explicitly in ( |6.1| ) that 
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M[l]uu 7^ if and only if \i G Pq.) The matrix if) is actually the submatrix of the usual 
S-matrix, restricted to Pq, and then rescaled by y/n. More specifically, it is given by 



Vv = c ex P 

where 



(n + k)n 



det 



y +[y] 

exp -2m — 

n + k 



(6.2) 



, — rt + 1 n(n.— 1) 

(n + /c)^i^^, (6.3) 



and as before t + (fx) = n(n 2 ~ 1 ' > + ^ i/x^, and = n — z + Yl^Zi Pj- 

For any A G P+ and \i,v G Pq, this NIM-rep coefficient J\f% equals the fusion coef- 
ficient Nji^ , where i is the unique solution (mod n) to the congruence ki + £(A) = 0. 
To see this, note that this i is the unique number for which J l \ G Pq. The ^-matrix 
factorises into the tensor product of the matrix in ( |6.2j ), and a U(l) level n S'-matrix of 
phases. Now, whenever the S'-matrix decomposes into a tensor product, then the fusions 
decompose into a product. So the su(n) level k fusion A x \i here will be J~ l ((J l X) x /x), 
where the fusion (J 1 A) x fx of weights in Po can be obtained by putting ip directly into 
Verlinde's formula ( |2.2| ). This means that ( |2.12| ) here will equal Nj ix ^, as desired. 

Hence we have proven that this choice of ip will always define a (non-negative integer) 
NIM-rep with the correct exponents. A similar analysis applies for instance to so(2n + 1), 
so(2n) and sp(4n + 2), when the level k is odd. 

6.2. A case with fixed points 

Next, let us consider su(p) with p > 2 prime, when p divides the level k (for p = 2 we 
require 4 to divide k). As we shall see shortly, that the formula we are about to give is 
indeed a NIM-rep follows from the conjectured formula in Section 6 of [^2| for the S'-matrix 



S e of the chiral extension by a simple current, as is also for example mentioned in ||11|| . 
As long as p 2 does not divide k + p, we can actually prove that our ifj matrix yields a 
(non-negative integer) NIM-rep. 

The exponents are the weights fx = (/j,q; fix, . . . , fi n -i) with p dividing t(p) = J2i Wi- 
There is a single fixed point <fi = jjj, . . . , , which has multiplicity p (as an exponent). 
The 'boundary states' consist of all J-orbits [u]. Again, the fixed point has multiplicity 
p. Although it may not be completely obvious, these two sets always have the same 
cardinality, namely 

1 Hn + p- 1)1 , \ . . 

' ' +p 2 -l ) . (6.4) 



p \ (n — l)\p\ 
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Define Vv,M> wnen neither \i nor v are fixed points by 



(6.5) 



where S 1 is the S'-matrix of su(p) level k. Of course this formula is independent of which 
representative v is taken for the J-orbit [v\. 
When n is the fixed point <p but v is not, then 



(6.6) 



for any z = 0, . . . ,p — 1. 

When fi is not the fixed point, but v = (f), then 
for any j. 

Finally, when both fx = v = <fr, we get for p odd 

1 



5. 



^,h),m,j) = - (x-i+ P 5 hij ) 



where 



x 



P 



k + p 



p-i 

2 



s for the Legendre symbol s 



(k +p)/p 



P 



(6.7) 



(6.8) 



(6.9) 



Recall that the Legendre symbol ^ ( fc +p)/p j equals if p 2 divides k+p; otherwise it equals 
±1 if (k +p)/p is or is not a perfect square mod p, respectively. So (|) = 1,0,-1 for 
a = 1, 0, 2 (mod 3), respectively, whereas (|) = 1, 0, —1 for a = ±1, a = 0, a = ±2 (mod 
5), respectively. 

On the other hand, for p = 2, the formula becomes 



-14 



k + 2 



(6.10) 



In this special case, the formula agrees with Eq. (B.6) of ||. 

In order to calculate with and analyse this ip matrix, it is convenient to use 'fixed-point 



factorisation' [33]. This is a way of computing S'-matrix entries involving simple current 
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fixed points. For instance the entries will equal either or ± [j^, J ■ The quantity 
x in (|6.9| ) equals S^. From this we quickly get that our will indeed be unitary. 

It is now straightforward to calculate the NIM-rep coefficients N\ x in terms of the 
WZW fusion rules. We find (whenever fi, v are not fixed points) 

V 

(6.11) 





= E JV - 

J 


yV A [<f>,i] 






= 1(n 

p V ' 



+ (p Sij - 1 



S\(j, 



\<j> ' W U IJ X I Q 



00 

Thus the first two are manifestly non-negative integers. We know (from fixed-point fac- 
torisation) that the ratio here will be or ±1. So it suffices to verify, for all A G P+, 
that 

NU = ^ ( mod P) • ( 6 - 12 ) 

>->O0 

Condition ( |6.12| ) is certainly true when p does not divide t(\) - - in that case we 
have Nf^ = = |r^-. When p does divide t(X), then our precisely equals the 

extended fusion coefficient ^[a]^]' f° r the conjectured ^-matrix S e for the simple current 
chiral extension, as given in |32|. This observation is also made in [11]]. So provided their 
conjecture is correct (as is expected), then the ip given above does indeed define a (non- 
negative integer) NIM-rep. Some relation between NIM-reps and extended S'-matrices is 
expected — see e.g. section 6 of []23|] , Thm. 4.16 of [ |T6"|| , and (B.6) in 0. 



S o<p 



Assume now that p 2 does not divide p + k. Then p does not divide the integer 
= + l\ . So it suffices to show that p must divide the integer 

1 ( M 4> S\$\ _^ s ^ s lf, S X( f> ( Sjcj, 1 \ , s 

^00 V ^00/ ^^O M ^O0 b o<f> \o ^ b 0(f) J 

where the sum is over all non-fixed-points \i. 
Consider the sum 

>o» s o4> [m] Sou 

where the first sum is over all non-fixed-points fj,, and the second sum is over all J-orbits of 

non-fixed-points \i with ^ 0. Note that this sum over will be an algebraic integer 

s 

(since each ratio is an algebraic integer, being the eigenvalue of an integer matrix) 
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V V ^ = pV^t (6 14) 



and it will be rational (since it will be fixed by all Galois automorphisms), and therefore 
it must be an integer. So p divides the first term on the right hand side of ( |6.13|) . 
Also, Fermat's Little Theorem implies that p must divide the integer 

sL 1 fk y- 1 , . 

1- - + 1 . (6.15) 



Thus p divides ( |6.13|) , and hence Af will indeed be a non- negative integer, provided only 
that p does not divide ^ + 1. 



7. Conclusion 

In this paper we have made a proposal for the boundary states of the WZW models 
that correspond to the conjugation modular invariant associated to a symmetry uj of the 
unextended Dynkin diagram. The boundary states of this theory are labelled by the un- 
twisted representations, and they are described by a formula ( |3.13|) similar to Cardy's 
expression for the boundary states of the diagonal theory. The corresponding NIM-rep 
agrees precisely with the twisted fusion rules. These fusion rules are given by a Verlinde- 
like formula (|3.14| ). 

We have given explicit formulae for the ift- matrix for all possible cases (section 4), 
and we have related the NIM-rep entries to untwisted fusion rules, thereby proving their 
integrality. We have also checked nonnegativity explicitly for numerous cases. Some of the 
corresponding fusion graphs have been collected in the appendix. 

Our formulae seem to generalise further to certain other classes of modular invariants, 
in particular those associated to simple currents. In order to illustrate this point, we have 
described two classes in detail in section 6. The results presented there already suffice 
to determine the NIM-reps for all simple current modular invariants for su(p) level /c, 
when p > 3 is prime. We have also proven that our NIM-rep coefficients are non-negative 
integers (at least when p 2 does not divide k + p) . The general case (that will be described 
elsewhere) will essentially cover all remaining modular invariant WZW models. 
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Appendix A. Explicit descriptions of the fusion graphs 

In this appendix we give explicit descriptions of some of the NIM-reps we have found 
in terms of the corresponding graphs. For the examples under consideration, the NIM- 
reps are often characterised in terms of the matrix associated to the field A = Aj. The 
corresponding graph has M vertices labelled by the rows (or columns) of J\f\, and the 
vertex associated to i and j are linked by (A/a) ij lines. Incidentally, given our discussion 
of section 3 above, these graphs are precisely the fusion graphs that describe the fusion of 
the field A with the twisted representations (that label the rows and columns of A/a)- In 
particular, the NIM-rep matrices considered here are therefore symmetric, and the relevant 
graphs are not oriented. 

For the case of su(n) with uj being charge conjugation, the graphs are explicitly given 
in Figure 1. The left- right symmetry of the su(4) graphs is due to the symmetry J B of the 
A y z> = D ( z> Dynkin dia gram. In addition, the graphs for su(4) are 2-colourable because 
of the symmetry J E of the diagram. Finally, there is also a left-right symmetry of 
the su(3) and su(5) graphs when the level is odd, due to the = A± and Dynkin 
symmetries J B , respectively. 

For the case of so(2A) with chirality flip, the relevant graphs are explicitly given in Figure 2. 
The graphs for so(6) do not agree with those of su(4) given in Figure 1 since, under the 
identification su(4) = so(6), Ai of su(4) corresponds to A3 of so(6). All of these graphs 
have an order-2 symmetry, which is due to the symmetry J B of the and Dynkin 
diagrams. The graphs are disconnected because of the symmetry J E of the and 
A% ] dia grams. In particular Q E (Ai) = 0, and thus the grading (|2.16|) implies that each 
component consists of the a G B with a fixed value of the grade Q E (a). 
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Fig. 1: The NIM-rep graphs for su(iV) with charge conjugation for N 
and k = 1, 2, 3. 
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2 2 2w 
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V 2—2 2^^-2 2—2 

Fig. 2: The NIM-rep graphs for so(2iV) with chirality flip for N = 3, 4, k = 1, 2, 3. 

For the case of so (8) with triality, and the case of Eq with charge conjugation, the relevant 
graphs are explicitly given in Figure 3 and Figure 4. respectively. 



2 21 





Fig. 4: The NIM-rep graphs for Eq with charge conjugation for k = 2,3, 4. 
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Fig. 3: The NIM-rep graphs for so(8) with triality for k = 2, 3, 4, 5, 6, 7. 
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